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INGHVI1 


IS 

o   principal  objective  in  the  design  of  a  blast-resistant 
structure  is  to  protect  the  structure  itself  including  its  e- 

ipment,  contents,   1   occupants  from  the  various  effects  of 
atomic  weapons •  The  design  consists  of  the  deto     tion  of 
the  load  acting  on  the  structure  as  a  function  of  tine, 
structural  resistance  re  quired  to  ms  of  the  in- 

dividual members  within  a  prescribed  maximum  value  is  calcu- 
lated by  dynamic  analysis  of  the  system.    0  limits  of  allow- 
able deflections  depend  on  economical  factors. 

A  numerical  example  presents  the  design  of  the  important 
elements  of  a  winclowless*  one-story,  reinforced  concrete  frame 
bull  i  •       last  loads  on  the  frame  are  calculated  as  t    « 

sted  by  the  Uni     bates  Atomic  Energy  Commission  and  the 
United  States  Army  Corps  of  engineers.  3 reliainary  design  of 

bers  i3  done  using  an  idealized  straight  line  load-time 
curve  and  is  checked  by  numerical  integration  using  the  calcu- 

-ed  load-time  data* 


■v    lu&te  student,      tment  of  Civil  *&:  .  . 
Kans<-     .   niversity,  Manhattan ,  Kansas. 


ISTaODUCIIOK 

The   three  most  important  phenomena  associated  with  an 
above-ground  explosion  of  an  atomic  bomb  are  air  blast*  ther- 
mal radiation,  and  nuclear  radiation.   In  designing  protective 

construction*  the  dynamic  load.,     ,  ed  by  air  blast  pressure 

2 
is  of  primary  concern  9  though  the  radiation  phenomena  must  be 

considered  in  the  design  of  structures  &     ed  to  offer  protec- 
tion for  personnel. 

It  is  not  possible  to  protect^  a  surface  structure  from  a 
direct  hit  of  any  size  of  atomic  or  hydrogen  bomb;  but,  there 
are  l^rge  areas  surrounding  an  explosion  in  which  conventional 
structures  would  collapse  or  suffer  severe  damage  while  a  blast 
resistant  structure  would  suffer  lamage*   be  cost 

of  construction  that  will  provide  this  protection  will  depend 
on  the  type  of  bomb  assumed,  the  assumed  disterce  from  the 
ground-zero  Cthe  point  on  the  ground  directly  below  the  explod- 
ing bomb),  and  the  degree   of  damage  that  can  be  tolerated  in 

4 

the  structure  • 

3 

/oral  types  of  construction  may  be  used  for  protection'. 
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esumably,  the  greatest  level  of  protection  can  be  secured  in 
buried  or  serai-buried  structures.    r»n  it  is  not     ^ble  to 
construct  a  buried  or  semi-buried  structure,  one  of  the  several 
forms  of  a  surface  structure  may  bo  considered,  such  as  shear 
wall  construction,  arch  and  dome  construction,  or  rigid  I 
construction.    adowless  construction  is  preferable  for  blast- 
resistant  structures,  because  it  eli     <js  the  danger  of  per- 
sonal injuries  from  nuclear  and  thermal  r.     ion  and  from  bro- 
ken glass. 

In  the  design  of  a  structure  capable  of  resisting  blast 
forces  which  are  large  in  magnitude  and  dynamic  in  character, 
members  and  joints  are  permitted  to  deflect  and  strain  much 
further  than  is  allowed  for  usual  static  loads.   If  the  struc- 
ture is  designed  to  resist  the  dynamic  blast  forces  with  stress- 
es in  all  structural  members  remaining  within  the  elastic  range, 
the  resulting  structure  will  be  more  expensive  than  a  structure 

in  which  plastic  yielding  of  a  reasonable  amount  is  permitted. 

x 
ffhe  amount  of  plastic  distortion  permitted*'  must  be  kept  small 

enough  to  provide  a  margin  of  safety    Last  collapse  and  to 
limit  damage  of  building  services. 

The  dynamic  character  of  loading  coupled  with  the  plastic 
yielding  of  members  re  Auires  that  the  design  procedure  be  based 
on  dynamic  anal/sis. -^    ier  the    ,1  rates  of  strain  that  oc- 
cur during  blast  loading,  materials  develop  higher  yield  stresses 

1  ll""*'— ^— w*— *^— im»mmm*mmm~~W9mmmmm***mmmm»mmmm*mm+*mmm~i^~Mm**»~*mm*mm*mm*KmMw*m<i'^  wmw»whww>wi>w»i  muni  nwi  hi— nw 
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orris,   op.    cifc.,    p«   2^6. 
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Army     c         of  engineer*,      •      .    Government  Prlntij 
March,   1957. 


than  they  do  in  the  case  of  statically  loaded  members.   The  in- 
creased dy   La  yield  stresses  are  used  to  determine  plastic 
strength  of  dynamically  loaded  members, 

WSAF0S~£FFSQ2  DATA 

explosion*  in  general f  results  from  the  very  rapid  re- 
lease of  a  large  amount  of  energy  within  a  limited  space.   he 

i 
sudden  release  of  energy  causes  a  considerable  increase  of 

temperature  and  pressure,  so  that  all  the  materials  present  are 

converted  into  hot  compressed  gases.   ince  these  gases  are  at 

very  high  tei  ,ures  ar    assures*  they  expand  rapidly  and 

thus  initiate  a  pressure  wave  called  a  shock  or  blast  wave. 

st    o  in  an  Infi  .         ecus  fcaospnere 

It  is  necessary  to  evaluate  various  aspects  of  air-blast 
phenomena  s  Lated  with  detonation  in  air  of  an  atomic  bomb 
in  order  to  determine  air  blast  loading  on  a  structure. 

Im     jely  after  the  occurrence  of  the  detonation*  hot 
gases  initiate  a  pressure  wave'  in  the  surrounding  air.    Lg* 
1  (a)  shows  '     ?neral  nature  of  v    I'lon-**'   of  the  air  over- 
pressure, i.e.,  excess  pressure  over  the  atmospheric  pressure. 
As  the  pressure  wave  is  propa     :  away  from  the  center  of  the 
,  the  following  (or  inner)  part  moves  th    h   a  region 

of  Mate'  -'uctural  Elements '' ,   .   • 

Arm.    vs  of     risers,   .   .   ov't.  J  rintii    r-'ice,  ftar.,1957» 

'  ffect  Data",  ,,    .  Lneers,  U«  . 

Gov      it   'inting  /'lice,  July,  1959*  paragraph  5. 


which  has  bean  previously  compressed  and  heated  by  the  leading 
(or  outer)  |  . ts  of  the  wave.  j  to  the  increase  in  the  tem- 
perature and  pressure  of  the  air  through  which  the  wave  is  rov- 
ing, the  inner  part  of  the  wave  iaoves  more  rapidly  and  catches 
up  with  the  outer  part  as  shown  in  .  1  (fc.  .  Thus,  the  wave 
front  gets  steeper  and  steeper  and  with..    very  short  period 

oe comes  abrupt  as    .n  in  Fig.  1  (c).  At  the  advancing 
front  of  the  wave,  there  is  a  very  sudden  increase  of  pressure 
from  normal  atmospheric  to  the  peak  shock  pressu      e  shock 
front  thus  behaves  like  a  moving  wall  of  essed  air. 

As  the  blast-wave  travels  in  the  air  away  from  its  source, 
the  overpressure  at  the  front  steadily  decreases,  and  pressure 
behind  the  front  fails  off  in  a  regular  manner.    ter  a  short 
tlmef  when  the  shock  front  has  travelled  a  certain  distance  from 
the  fireball,  the  pressure  behind  the  front  drops  below  that  of 
the  surrounding  at    iere  forming  the  negative    se  of  blast- 
wave  as  shown  in  fig*  2.   The  symbol  ;   denotes  the  peak  over- 
pressure in  pounds  per     Pi     ,   nd  U  is  the  velocity  of 

shock  front  in  feet  per  second. 

o 

Figure  3  shows'  the  varietion  of  overpressure  versus  time 

of  trave  .  The  symbol  t  denotes  the  duration  of  the  positive 
phase.  •  overpressure,  i_,  at  any  time  t  after  the  arrival 
of  the  shock  front  is  given  by^f '  the  expression: 
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FIG.  1.   OVERPRESSURE  DISTRIBUTION  IN  EARLY  STAGES  OF 
SHOCK  FORHA.TION. 
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caling  the  Blast  i henomena 

7 
It  has  been  found-'"  that  air-blast  phenomena,  such  as  the 

pressure  and  duration  at  different  distances,  are  related  for 
different  strength  boobs  according  to  the  ratio  of  the  cube  root 
of  the  equivalent  weights  of  TNT.   These  relationships  are  re- 
ferred to  as  the  scaling  laws* 

These  scaling  laws  stete  that  if  a  giver  peak  overpressure 
is  experienced  at  a  distance  r,  from  an  explosion  of  bomb  of 
total  energy  yield  tf,,  the  same  peak  overpressure  will  be  expe- 
rienced at  a  rg  from  the  explosion  of  a  bomb  of  total  energy 
yield  W*2,  where'' 

~  -  (~ )  (2) 

The  same  scaling  laws  also  state  that  while  the  peak  pres- 
sure from  the  two  bombs  are  equal  at  the  two  radii  r,  sad  r~, 

the  durations  of  the  blast  pressure  wave  at  t  Q  oints  are  dif- 
ferent. If,  for  example,  the  duration  of  the  positive  phase  of 
the  pressure  wave  from  the  first  bomb  is  t  ,  at  a  distance  r^, 

the  duration  of  che  positive  phase  of  pressure  wave  from  the 

7 
second  bomb  at  distance  r~  will  be', 

Loading  on  structures 

The  manner'  in  which  the  blast  wave  loads  a  structure  is  a 
ction  of  distance  of  the  structure  from  ground  aero, 
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FIG.  4.  A  BL, 

Bt  A  STRUCTURE. 


Under  the  condition  that  the  blast  wave  has  not  yet  com- 
pletely surrounded  the  structure ,  there  will  be  a  considerable 
pressure  difference  between  the  front  and  back  face...    jh  a 
pressure  difference  will  produce  a  lateral  force,  tending  to 

ise  the     'ture  to  move  bodily  in  the  same  direction  as  the 
blast  wave.   This  force  is  known  as  the  "difl     ;n  loading". 
b.        oadinff.   rag  loading  is  the  term  given  to  the 
forces  on  a  structure  resulting  frois  the  fa      Locity  of  the 
air  :.ar~icles  in  the  air  blast  acting,  as  a  hi     leci : 

pe  of  loading  is  most  important  on  truss-ty^e  structures 
such  as  bridges,  and  buildings  in  which  the  wall  panels  fail, 
leaving  the  structural  frame  exposed  to  the  air  blast* 


ALY3IS  FOR  DYNAMIC  LOADS 

Q 

rwo  different  methods  are  used     either  separately  or 

8 Principles  of   Jynamic  Analysis  and  JJesicn",      •    •      w 
ngineers,   U«    .      ov*t.   printing  offioe,   inarch,   i;;>7» 


10 


eoncurr  er  d, 

on  t  ik 

do:         »rgy  c  ii- 

librium  of  Kewtoi        Ion  of       , 

>n : 

I    -----  a)  (4-) 

to  the  prii      of  cor;  , 

nation  i&  ot 

>ne  o     :ic  energy  (5) 

in  t    Lon  (;,  rev< 

in  oner  ;,n 

iitt  :€  un- 

der lc  j.xje, 

i.e. 9  (aaa 

i  in  the      ce  of  kinetic  e 
energy  cone.    Lon«  to  the 

•    nee,  a 

^t   co  in  th 

fy  3tea 

In  d„  be 

coi:  : >  the  strain  energy,  ic 

;  ay 

int  c 

is  c  in 

i  the  structural  eJ      ,  ,  .  • 

el   ,  on 


XI 


and  rot    n  of  all  the  masses  of  the  structure.   The  actual 
evaluation  of  these  quantities  for  a  given  structure  under  dy- 
namic load  is  very  complex*   For  this  reason,  it  is  necessary 
in  many  problems  t  i  11m  '-oth  the  structure  and  the  load- 

ing,   a  distributed  masses  of  the  given  structure  are  lumped 
together  into  a  number  of  concentrated  masses,  and  the  strain 

rgy  is  assumed  to  be  stored  in  a  number  of  weightless 
spri) 

when  idealizing  dynamic  load  on  structures,  I  p   impllfi* 
e  generally      red.    is  first  one  involves  the 
geometric  distribution  of  the  load  over  the  structu:  «   If  for 
the  purpose  of  analysis  the  mass  of  the  system  is  concentrated 
at  certain  points,  the  lo    uat  be  applied  at  the  same  points. 
This  requires  nodi     ion  of  the  magnitude  of  the  lc 
second  simplification  involves  idealising  the  lc  .- ■  ,-ime  curve. 
If  a  numeri     ethod  of  analysis  is  used,  it  is  not  necessary 
to  idealize  this  function  I       any  variation  oai  >•  handled. 

.■•.-,  ;ts  merely  of  a  number 

of  concent  Lned  together  by  wei-^htless  i     s 

I  subjec  aitrated  loads  which  vary  with  time,   fig- 

ure 5  shows  two  simple  structures  together  with  corresponding 
idealized  or  equivalent  dynamic  systems. 
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"  •  6  i:  8,^\  in 

1  Lch  r  (t)  ±enot<  I  ic         ttloa  of  timo, 

m^,  t/    3s  of  ei'        ,  *#Wi    Latomal  resia- 
»•  fl  '.x  is  the  defl  I   km* 

Application  of    ton9!  o  oatloo  of  notion  .fields  t 

d2x 

where  — *  denotes  the  accD_    Ion 

dtr 

Equation  (?)  -jives I 

"  Pe(t)  dx  =  I  ne   (ff)2  ♦  $  X  B  Oc)  dx        (7) 
o  o 

dx 

in  which  gx  denotes  velocity 

e  point  of  maximum  deflection,  i.et,  ze^o  velocity, 
Lon  (7)  yields: 

{  **   P  (t)  dx  -  yXm  R  (x)  dx  (8) 

o  o 

where  x„  ■  maximum  deflection 

us 

1  ork  Donoi 

A  typical  dynamic  lo:  4  i  i   shown  in  Fig.  7»  where  I  is 

defined  as  the  duration  of  the  lo  ..d.   If  this  load  is  applied 

to  stem  shown  in  Fi~.  6,  W  (t),  the  external  work 

e 

done  up  to  any  tiao  t,  is  ;\iven  byt 

r   x(t) 
VL(t)  -       PAt)   dx 


£'  (t)  ox.  dt  (9) 


e      '0 


o  *•"'» 
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Integration  of  equation  (6)  yields: 


&  .  1 
at   m 


e  J  o 


Qp6(t)  -  &#oo 


at 


(10) 


?husf  the  expression  for  the  work  done  becomes 


v*> 


p.ct) 


(  *•<*)  -  HL(x))  ** 


e  o 


dt 


P(t) 


FIG.  7.   PYiTCAL  JOTAHIC  LOAD 


At  time  t  corresponding  to  the  maxiroia  displacement  of 


the  mass,  the  maximum  work  done,  W f  is  given  by: 


TMI 


■ 


0 


p.ct) 


e  ^o 


(  Pe(t)  -  Re(x))  dt 


dt   (11) 


If  t  is  much  greater  then  3?.  the  resistance  H  (x)  is 
m  °  '  e 

small  during  the  application  of  load  and  may  he  no  -  leetod. 
Thus,  the  expression  for  work  done  becomes: 


*9* 


f 


Vt} 


J-  [  t   *e(t)  dt 


e  o 


ut 


Where  W_  denotes  the  work  done  ignoring  the  contribution 
of  the  resistance. 
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Integration  of  equation  (12)  yields: 
in  which 

H^  is  the  tot, 1  impulse  of  the  external  load  and  is  equal 
to  the  area  under  the  load-time  curv  . 

f ter  the        - m   of  the  load,  the  mass  has  acquired  a 
kinetic  energy  equal  to  the  work:  done: 

H2 


srS  -sme  ** 


Therefore,  the  initial  velocity,  Vf  is  given  by 

v  -  %  (W 

es,  the  internal  resisting  force,  ?L(x),  cannot 

e 

be  neglected  in  the  time  interval  between  0  and  T#  Therefore, 
the  work  given  by  equation  (13)  may  be  considered  to  be  the  ab- 
solute maximum  work  which  could  be  done  by  a  given  load  on  a 

W 
dynamic  systen.   i'he  ratio,  rjF^  ,  of  the  actual  work  done  di- 
ps 
vided  by  the  maximum  work  is  called  the  work-done  ratio. 

Curves  giving  the  work-done  ratio,  Cw,  for  two  simple  load- 

time  functions  and  for  one-degree  of  freedom  are  plotted  for 

design  purposes.   Theas  are  given  in  terms  of  ratios 

8Ibid.,  Fig,  5*24  to  5*2?,  p.  50. 
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°R  *  "IP  "■*  GT  "  T"  (14a) 


where  8^  denotes  the  maximum  or  plastic  resistance  of  the  sys- 
tem, B  the  peak  value  of  external  loads,  T   the  duration  of 


I  M 

f5-  is  the  natural  period  of  the  system* 


Transformation  Factors 


The  application  of  transformation  factors  to  the  dynamic 
parameters  of  a  structure  transform  the  system  to  an  i    Lsed 
system.^  they  are: 

(a)  Load  Factor  1U.  The  concentrated  dynamic  load  on  the 

equivalent  system  is  obtained  by  multiplying  the  total  load  on 

the  actual  structure  by  the  load  factor. 
P_(t) 

H  -  -Ht7  (15) 

K»  is  determined  by  equating  the  external  work  done  by  P 
on  the  equivalent  system  to  that  done  by  the  real  syster  . 

(b)  flees  Factor  Ku..   When  the  total  mass  of  the  structure 

is  iLultiplied  by  the  mass  factor,  the  concentrated  mass  of  e- 
ivalent  system  is  obtained. 

hi  -  s;  <i6> 

i'his  factor  is  obtained  by  equating  kinetic  energy  of  the 
real  system  and  the  equivalent  system. 

<uorri3,  op.  cit.,  p«  149. 
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(c)  Resistance  Factor  K    The  resistance  of  an  element  is 

the  internal  force  tending  to  restore  the  element  to  its  equili« 
brium  position.   At  a  given  deflection,  the  resistance  is  de- 
fined as  being  numerically  e^ual  to  the  static  load  required  to 
produce  the  same  deflection,   ihe  product  of  the  resistance 
factor,  K~,  and  the  resistance  of  the  real  element  gives  the 
resistance  of  the  spring  in  the  equivalent  system.  Equating 
the  internal  strain  energies  of  the  two  systems  yields: 

-  -§  (17) 

(d)  iiaximum  Resistance,  IL  and  Spring  Constant  H  .  The 

maximum  resistance  of  a  real  element  is  defined  as  the  maximum 
total  load  which  can  be  carried  by  this  element.    e  product 
of  the  res     ce  factor,  K«  and  the  maximum  resistance  give 
the  maximum  resist: nee  of  the  equivalent  syste.  • 

fhe  spring  constant,  R  of  the  real  system  is  defined  as 
the  total  static  load  to  cause  a  unit  deflection.     ice  the 
deflection  of  the  two  systems  should  be  the  same,  the  spring 
constant  of  the  equivalent  system  is  obtained  by  applying  the 
resistance  factor  thus, 

*me  '  KR  *m 

(18) 


(e)  Load  Mass  factor  K^.   inhe  ratio  of  mass  factor  to 
lot         l     ined  as  the  load  mass  factor. 


la 


(f)  ftynasxlc       pa  ?•  Dynaiaic  reactions  are  needed  to 

determine  the  shear  at  the  i  ^ts  and  are  obtained  by  v 

oration  of  the  resistance  of  the  e     lent  system  and  ©v 
applied  loads* 


nsfonaation  Factors  ported, 

Uniformly  loaded  Beam* 


As  an  m        '..e,  %t   .  :,ors  are  calculated  below 

for  a  sim;  1^*       ted  t    *         nifora  sass  and  subject- 

ly   distributed  load.    ig,  B») 


P(t)  »  p(t)  L 


t t i *  v  v i x i i J 

^    b^  ■  m  L     ^ 


M 


(b) 


e  e  e 


(c) 


FIG.  8. 


i  ■  J 


(c,        ST 


objected  to  ds% 

rps  of    -n;  in-.;-;    v      .      ,       weswM&t  Printing  Office, 
-   i    I   57 1  p. 
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(a)  Load  r actor*   The  elastic  curve  of  a  simply  supported 
beam  subjected  .    uniform  static  load  p(t)  lfcs/ft  is  given 
by  i 

xa  -  ^-j   (L5  -  2Lf*  7?)  (20) 

The  deflection  at  the  mid  span  is  given  by 
.4 


xac  "  3 


Thus, 

xa  *  *%  (I,3z  "  2L*3  *  ^  xac  (22) 

2he  totul  work  done  by  the  load  is  equal  to: 
*a  "  V.   »-«•  d4 


Id     ac 


(23) 


where  3?  is  the  total  load 

In  the  equivalent  system,  the  work  done  by  the  external 
load  is, 

u     e  e 

ing  x  equal  to  x   and  equating  the  work  done  JL:  ooth 
eases,  gives: 

*.  •  i " 

Thust  the  elastic  load  factor  is, 

After  the  formation  of  plastic  hinge  at  mid-apan,  the  de- 
flected shape  is  assumed  to  be  as  shown  in  Pig*  9* 
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P(T)  «  p(t)  X, 


}    1  j  1  1  1  II  J 


k-  Z 


Plastic  hinge 


v2— H 


(•) 


He-Bw 


J  *© 


(b) 


fig*  9. 

.   (a)  (i         OT 

- 


In  this  case  K,  »  4 


(25) 


(b)  fS&sp  factor*  The  maas  factor  is  obtained  toy  equating 
the  total  kinetic  energy  of   the  beam  to  that  of  an  equivalent 
eyste  ,   In  the  simple  harmonic  motion  of  the         •  assumed  heref 
the  maximum  velocity  at  any  point  along     beam  is  proportion- 
al to  the  ordinate  of  the  deflection  curve  at  the  same  point* 
,  in  the  elastic  range , 

(26) 


a   51* 


*e 


where  V   is  the  velocity  at  mid-  n   and  V  is  the  velocity  at 


M 


a  distance  Z  fro&  the  support. 

The  kinetic  energy  of  the  real  system  is  given  by: 
L  2 


I 


Va     d* 

a 


(27) 


where  m  e  mass  per  v,nit  I 


■1 


Integrating  equation  (27) •  yields: 

KEa  -  0.25  *i   Vac2  (28) 

where  m.  is  the  total  mass  of  the  beam* 

rhe  kinetic  energy  of  the  equivalent  system  is  iven  byj 

-  |  *«  ve2  (29) 

Equating  the  kinetic  energies  of  the  real  system  and  the 
equivalent  system,  yields  the  elastic  mass  factor. 

In  the  plastic  range,  the  deflected  shape  is  shown  in 
•  9  (a)  and  the  mass  factor  is  given  by, 
%  -  0.33  (3D 

.C;   flaxinma  Resistance  and  opring  Constant.    e  maximum 
moment  for  a     ly  supported  beam  at  mid-span  is  given  by 

-'■c  '   ^ 

where  i  is  the  sum  of  the  uniform  load.    I  maximum  resistance, 

H  is  obtained  by  equating  the  bending  moment  to  the  resisting 

moment.   Thus, 
8  fl 

B^  m   1£  2) 

In  the  tqul         3tem,  the  maximum  elastic    i3tance 
is  given  I 

This  is  the  lisiitl'   -ssistance  in  the  elastic  range*       ;r, 
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in  the  plastic  ran  <  ,  the  resistance  of  the  equivalent  system 
is, 


i   8  M 


(5*) 


l«  stiffness  or  the  spring  constant  of  the  real  beam  in 
the  elastic  range  is, 


k.MijH 


and  the  spring  constant  of  the  equivalent  system  iss 


(35) 


k  - 


v  16  m  m 


KLk  B  "5 


(56) 


(d)  Ifrnamic  e actions *  In  order  to  determine  the  dyna:.  ic 
reactions  on  the  actual  beamt  the  inertia  forces  distributed 
along  the  beam  must  be  considered*  In  Fig*  10 ,  it  is  assumed 
that  the  inertia  forces  are  at  all  points  proportional  to  the 
ordinates  of  the  deflected  shape  of  the  beam* 


A 


V 


P(t)  -  p(t)  L 


1  1  1  1  1  1  llB 


Kt) 


,   u  j  2£|i 


61  L  %   I(t) 

h-192  ~* 

h — H 


(a) 


Assumed  distribution 
of  inertia  forces 


(b) 


cs*  io.  i  ac 

HANGS,   (a)  LOAD  ABP  U      -FORC.  01  ■;  (b)        OK 

0K&-  . 
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Considering  one-h^lf  of  the  beam  and  taking  moments  about 
point  i)t  the  centroid  of  the  inertia!  forces,  gives: 


v^l 


m   1}c  -  \  p  (lf§  L  -  L/V  "  °  (*7) 


where  H  is  the  b<  Ling  moment  at  mid-span, 

8  K 
Assuming  the  resistance  is  equal  to  — *&   ,  and  substituting 

for  li  in  equation,  yields: 

V  -  0.59  SL   «►  |  P  (; 

Using  the  same  procedure  in  the      .c  rangs  along  with 
the  assumed  deflected  shape  for  this  range,  «?ig«  9  (a),  yields: 

v  »  |  R  f  i  P  (39) 


AKALI  IS  i3T  KU1- 

Ihe  analysis  of  a  dynamic  system  of  a  single  decree  of 
freedom  consists  of  the  evaluation  of  displacement  using  Kew- 
ton's  e  >uation  of  motion  (i.q.  #>)*  Uwo  approaches  i^re  general- 
ly used  in  solvin  this  type  of  differential  equation*    e 
first  one  is  the  exact  method  in  which  the  solution  is  obtained 
directly  from  the  differential  e  puution.   This  approsoh  is  ap- 

ble  only  whan  both  the  load  and  resistance  function  can  be 
wd  in  simple  i    -.mtical  forms,     e  other  approach  is 
the  numerical  int.     ion  which  is  chorally  applicable  to  any 
I  of  load     resistance  function, 

tmmm — — mmmi^m^^mm fmmm mm »   i   i      ■      ■ ■  i— *am— i^n       i  n    ■ won ■ I— — ^11   n  i  i  ■  n     ■■  m  I  ■ I ■      i  i  — «^ *mmmn mmf 

rinciplas  of  ilynamic  nalysis  and  Design",   ♦   .  rmy 
Corps  o.J    ineer  ,   .   .  -overnaant  Printil    ffloa4   POh, 

7. 
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dseiplea  of  Numerical  /dialysis 

The  nuaerical  met  od  of  evaluating  the  displacement  fro» 

£quatloa  4  is  called  the  method  of  numerical  integration*'* 

e  load,  resistance,  .acceleration,  velocity,  and  display* 

ments  are  plotted  versus  time  la    •  li  i  ••..  lI  q  e. 

oume  that  t  ,  t, ,  tol  ..••,  t  , ,  t_ ,  and  t„Ll  eonsti- 
o*  1  d*  *  n-l*  n'      n+1 

tute  a  ti*ue  sequence,  an    it   t  denotes  the  tin*  interv 
from  t  to  t  t»        wnio  load  la  ee    I  to  be  initiated 
at  t  *  t  •   The  acceleration,  velocity,         iace&ent  at  t 
are  denote,  bj  ^ ,  v  ,  end  x  ,  respectively.        nccelera- 
tion  in  the  time  interval  ^t  is  represented  by  a(t),  the  ve- 
locity and  dis:     ent  at  t  ,  are  f;;iven  by   the  following 

ions: 

t 
n+1  *  fn 


K^    •  ▼«  ♦  \   n+1  a(t)  dt  (40) 


ta 

and 


*n+l 


xn+l  -  xn  +  vn  (A  V  +  J        La(t>_  dt  -  dt     <41> 

These  equations  indicate  that  the  velocity  and  the  dis- 

t  n  can  be  ob;     1  by  extrapolation  from  the 

corresponding  values  at  tw,  once  the  acceleration  in  the  time 

n* 

interval  ^t     is  known. 


5 

,    Q    »cit*|  p.    13$,   chapter  8. 

o 

"Principles  oi  La  ,      .      .      rmy 

Corps  o  ineers,      .      .      overi:  rintii  rch, 

1957. 


f(t)    and 
R 


to       t.     t- 


th-  feh  fc*M 


(a)   Load  and  Resistance 
a 


(b)   Acceleration 


(c)    Velocity 


(d)    Displacement 
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f(t)-R 

m 


v  =  ^adt 


v  dt 


FIG.  11.   LOAD,  RESISTANCE,  ACCELERATION,  VELOCITY,  AND 
DISPLACEMENT  VERSUS  TIME. 


>6 


In  the  analysis  of  structures  under  dynamic  load,  the  ve- 
locity and  the  displacement  at  the  time  t  *  0  are  assuaed  to  be 
equal  to  zero,  that  is,  , 

V  «  0  and  x^  *   0 
o  o 

In  applying  ^iuations  (40 > and  (41 v ,  .he  values  jf  v,  and 

x,  can  be  obtained  provided  a(t)  is  Joiovn  in  the  ti    tterval 

from  t  to  t,  •         >cess  can  be     inued  until  the  v.  lues 
o     l 


of  v„  and  x„  for  any  value  of  n  are  obtaij 

n     n      ^ 


*5  ,< 


Many  extrapolation  formulas  derived^  for  i 

solution  of  differential  equation  by  assuming  a  ;     a  accelera- 
tion-time relation  i      tisie  interval  t  •     extrapolation 

x    ft 

formula  based  on  acceleration  impulse  -*•  is  described  below. 
.celerdtion  Impulse  uxtrapolation  Method 


n+i 


(a)  Actual  acceleration 

curve • 


o 

H 

e 

•«3{ 


I(tfi 


kVi) 


*a-l     tn     *»+! 


(b)     vCccle,  ft  pulses  to 

replace  the  p;iven 
I  ft  curve. 


^ris,   op.    cit.,    [...    186. 

;ibiu.,   ,„.   191. 


2? 


V 

— 

+ 

Vn 

-1 

v~ 

Ci_ 

n 

*n-l   &n   vn+l 
(cy   Velocity 


*n-l  **       *n*l 
(d)  Dlaplao«m«ft$ 


FIO.  12.  AC   SRAfXOl  IMPULSE  EXTRA     ION  MBTH03. 

la  this  a©-    ?  the  actual  ace  ,;n  in 

Fig,  12(a)  is  ;   t  ced  by  a  tr  in  o?  t  iuallj     .ad  impulse* 
occ      at  tQ,  t-p  fc2*  .  .  .  tn.         it;  ~e  of      eel« 
©ration  impulse  at  t  i3  given  by: 


I(t  )  -  a  (^t) 
N  xr    n  x   ' 


where 


A(t)  -  tn  -  t  *  t0  -  t-, 
1    o    2    1 


(42) 


%a  •  t«  »  «  .  »  *  t.rt  -  n-1 
p     <£  a 


.iince  an  impulse  is  applied  at  t  t  there  is  discontinuity 
in  the  value  of  velocity  at  t  .   In  the  time  interval  from  t 
to  t  , ,  the  velocity  is  constant  and  displacement  varies  lin- 
%Iy  with  time.   t~  and  t*  denote  the  time  S        /  before 

lie  !  LOB  of  the  impulse  at  t  f  and  /"  and  V* 
denote  the  corresponding  velocities  which  are  related  by: 


n    n  v   ' 


(*« 


The  relation  between  x^  ,  and  x  .  and  between  x„    l  SL... 

n— l      n  n     n+± 

are  ^iven  by: 
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Lona   (4^)   and   (44)  yieli  r- 

re-  a  &  formula  f 01  ion  i; 

,2 


x 


n+1 


2  **  "  Xn-I  •  an  <  At>  5> 


je  the  v  1   of  :  *n-l       Y'   t?  rnf   tlTe  value 


of  x  at  t     - 


n  exa         1   t  ;.   .    i 

to  Bt  of  ■  • 

i  zero. 
e 
at  motion. 

Xh 
and  in  .  low- 

Q 

: 

t*       m  oi 

c 

Bg  concrete  1  f*ac  s   1 

lulus  of  elasticity  of  concrete,   E_  *  5(1. 

Hatio  c  ..of  el          sity                                          of 

elasticity  of   .  srete,  n  «  10 

Q 

J  9  •  • 

y  Cor  g   oi  , 

Jul yf    1959. 
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reinforcing  steel  J 

uatic  yield  point  stress,  f  «  40,000  psi 

namic  yield  strength  of  steel,  f ,  -  £2,000  psi 

Load  Determination 

The  computations  for  the  various  presoure-tiiae  curves  are 
based  on  the  formulas  given  by  the  Corps  of  ^n^ineers,  reference 
7. 


The  peak  overpressure,  *aot  is  obtained  from  Fig.  1,  Appen- 


dix B 


pso  "  10  Psi 

The  duration  of  positive  phase,  t  ,  is  obtained  from  Fig, 


2,     idix  B, 

t  *  0,71  sees 

Caj  ,  Determination  of  Front  Face  Overpressure  Versus  fim* 
Curve. 

velocity  of  eound,  B,^,  in  the  re  .ion  of  the  refl.ot.4 
overpressure  is  obtained  from  Reference  2,  ig.  11  -  '1, 

°ref  1  *  12yc  ** 

Jhe  time  re  paired  to  cle;jr  the  front    .  of  the  structure 
is  Kiven  by  Equation  (11.4)  Reference  ;, 
6   h' 


c   Grefl 

*  0.0^49  sees. 
For  P       »  10  psi,   reflected  ov«.  re,  f ,,   is  given 

by  J?ig.   11.20,    reference  3» 
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*refl  •  25.3  psi 

i?or  xc,n   "  10  psif  dynamic  pressure,  q  is  -;iven  by, 

q  »  2.2$     si 

rage  overpressure  on  the  front  face  1    /en  by  the 
.tion 

i 


front 


»  Ps  ♦  0.35  q 


tine  • 


Table  1.   e  termination  of  front  fytee  overpressure  versus 


t 
sees 


t 


From  Fig.  3 
Appendix  B) 


i 

psi 


s 


r. 


so 

•  4 
Appendix  B) 


s 
psi 


front  - 
Pg+0.85q 

psi 


0 

,.1 
0.2 

0.3 
0.4 

0*5 

0.6 

0.71 


0 
.141 
0.282 

0.423 
0.564 

0.705 
.346 


0.525 

0.267 

.118 

0.060 

0.025 

0.008 


2.230 
1.170 

0.595 
0.264 
^.134 
0.056 
.  18 


1. 
J.  746 

0.543 

0.360 
0.247 
0.146 
0.066 
0.005 


1J.00 
7.46 
5.43 
3.60 

, 
1.46 

.  ■ 

.05 

J 


1  .12 
8.46 
5. 

3-83 

2.58 

.68 

.05 

J. 00 


+ 


Fig.  17  shows  front  face  overpressure  versus  time  curve. 


(b)   vvera  ;e  Back  r.all  tveri/resoure  Versus  ^ii^e. 
For  1        •  10  psi,  the  velocity  of  incident  shock  front,  U( 
is  c;iven  by  Fig*  11.21,  nefereixe  3t 
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U0  -  1403  psi 


The  time  displacement  factor,  t,  ,  is  given  byi 

L 
A  "ft 


t,  .  fc-  (46) 


o 

where,  L  is  the  length  of  the  structure, 
Equation  (46)  yields: 
tj»  *  0*024  see. 

Tine  required  for  overpressure  on  the  re&r  face  to  rise 
from  zero  to  its  maximum  value  is  given  by:*^ 

o 

where,  C  is  the  velocity  of  sound  in  air  and  is  equal  to  1115 
fps. 

Apptadia  B,  Fig,  4,  for  t  -  td  m   0*0?^6,  the  value 


s  is 

found  to  be j 

of 

so 

Jl 

so 

*  J ,859 

or 

s 

•  8,59  psi 

The  peak  vnlue  of  average  overpressure  on  t     ck  wall  is 
given  by  £qpatlo$  (11 »9)  Reference  3  s 


(pback>  ««-i8b  X  *   a*^>  •  (48) 


where 


sb  w  incident  blast  wave  overpressure  on  back  face  t 
time  t  -  t_  •  tv. 
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F- 


°»5  K„ 
so 

14.  r 


(i;t.  .  O  max  •  peak  value  of  average  overpressure  on  back 
wall  which  occurs  at  time  t  ■  t.  +  t^ 
i^uation  (48)  yields: 

(Pback)  *nax  •  6.30  psi 

For  tixse  in  excess  of  t  •  t,  •  t^»  overpressure  on  back 
wall  is  given  by      ion  (11.10    I arance  3. 


back 


ack)maX 


1  + 


(Pback)aaX~ 


f 


•% 


'*-«VV 


-  t, 


(49) 


where 


t  is  the  duration  of  positive  phase. 


Table  2. 


Determination  of  the  back  wall  overpressure 
versus  time. 


t 

t-td 

t-td 

5i 

*8 

*back 

Pback 

sees. 

sees. 

*o 

so 

(Fig.  4 
Appendix  B) 

psi 

h 

psi 

0.078 

0.054 

0.076 

0.859 

8.59 

0.734 

6.30 

0.095 

0.U71 

0.10 

0.814 

8.14 

0.734 

5.96 

0.166 

0.142 

0.20 

0.655 

6.55 

0.739 

4.84 

0.308 

0.284 

0.40 

0.402 

4.02 

0.766 

3.08 

0.450 

0.426 

0.60 

0.220 

.20 

0.820 

1.80 

0.598 

•?68 

0.80 

0.090 

0.90 

0.901 

0.81 

0.734 

0.710 

1.00 

0.00 

.... 

1.00 

0 

'A 


Sample  calculation  of  Pable  2. 
i'or      t  -  t. 


*   *  0.40 


t  •   (0.40)   tQ  ♦  td 

•  0.508  sec. 
Substituting  the  values  of  O^cJp  »ax  »  6.30,  Pgb  •  8.59, 

t  •  0*306,  td  ■     ,        ,       d  tb  *  0.0538  in  Equation  (49)  yields? 

— *?ack>    *  0.7C6 

or 

*back  *  766  (4.02) 
«  3.08  psi 
Lg.  18  shows  the  back  overpressure  versus  ti.     rve. 

(c)  *.et  Lateral  Overpressure  Versus  Ij&e.   it  any  time,  t 
ral  overpressure  is  given  by: 

Pnet  *  pfront  ~  Pback  (5°^ 

Figure  19    v;s  net  lateral  overpressure  versus  tisse  curve. 

(d)  Average  ..oof  Overpressure  Versus  Time.  Figure  15  shows 
the  relation  between  the  overpressure  ratio,  ^roof   au<i  the  time 

t  in  sees.' 


I   -H 


»«  t 


FIG.  15.        oof  c  Tins 

FGH  A  jiGUL*         IE. 


re 


.   ♦  C.l  U  -  if**) 


SO 


(51) 


Pn   should  always  be  smaller  than  1. 

i'he  valae  of  P*  should  be  taken  by  the  equations  (52), 

whichever  gives  a  smaller  value* 

P        i  1/3 
P'  -  2  -  (jfSy  ♦  !)(£) 


or 


(52) 


P1  «  0.5  ♦  -.125  (2-  -jf^)2 
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Eat  time  displacement  factor,  td  •  ^tr- 

o 


0.012  sec. 


At  time  t  »  t ,  f  -gsr- 


0.024  sec,  Equation  (51)  yield* i 


56 


%^£  -  0.91 


At  time  t  -  & 


0.120,  Rations  (J?2)  yield: 


^2^  .  0#720 


or 


roof 


JF»  «  0.71? 
s 


therefore 


P 
pt  .  |°°^  •  0,7X7 


At  time  t  »  §^  +  ^  «  0.28  sec, 


o    o 


P 
p°°   »  1.0 
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Table  , .   Je termination  of  the  average  roof  overpressure 
versus  time. 


t 

t-td 

a 

■. 

-. 

P 

roof 

p 
roof 

sees. 

sees. 

*o 

'.0 

(  igure  4 
endix  B) 

psi 

psi 

0.024 

0.012 

0.01? 

0.966 

9.66 

.  ;io 

8.75 

0.083 

.71 

0.10 

0.814 

8.14 

0.791 

6.45 

.120 

0.108 

0.154 

0.726 

7.26 

0.717 

5.20 

0.156 

0.144 

. 

0.655 

6.55 

0.??9 

5.10 

0.225 

0.213 

0.30 

0.519 

5.1 

0.901 

4.68 

0.296 

0.284 

0.40 

0.402 

4.02 

1. 

4.02 

0.567 

0.355 

0.50 

0.303 

3.03 

1.00 

3.03 

0.438 

.  ,26 

.60 

0.220 

.20 

ItOO 

2.20 

0.509 

0.49? 

0.70 

0.149 

1.49 

1.00 

1.49 

0.580 

0.568 

.-30 

0.090 

0.90 

1.00 

0.90 

0.651 

0.639 

0*90 

0.041 

3,41 

1.00 

.41 

0.722 

0.710 

1.00 

0.0 

0 

1.00 

0 

U'ig.  20  shows  aver&ge  roof  ovei\  curve.) 
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Design  of  -all  Blab 


The  wall  is  dt     1  as  a  one-way  reinforced  concrete  slab 

spanning  from  fixed  support  at  foundation  to  a  pim-ed  support 

i  ■ 
at  the  roof  slab.  v  Q -e  slab  is  permitted  to  deform  in  plastic 

region  by  developing  plasti*  hinges  at  tve  foundation  and  near 

*ld-heignt«   Uhe  span  of  slab  is  equal  to  clear  height  of  wall. 


P(t) 


PIG.  21.  IGK  OF  WALL  SLAB. 

(a)  Deai;ra  Loading.  Design  loads  as  idealized  from  the 
commuted  loading  shown  ti  *    17  is  defined  by: 

Bt  leak  value  of  load  I     *cl  to: 

(??.?,-  W  (16.» 

1G00 

m   60  kip 


"single  otory     rame   Buildings".      .      .      rcay  Corps  of 
Lneers,   U.      .     oven  sent     rintine         ice,   Jan.   19u.    . 


*  •  *■ 


Duration  of  external  load,  T  *  0.062  sec. 
Pulse  Loading,  H  is  given  by,  II  »  B  lyp 

*  1*86  kip/sec. 

(b)  I'he  Avn.-^ic  i'-asiKn  .•actoiv*  e  dynamic  fac- 
tors for  the  slab  sh  1  •  21  are  obtained  from  Appendix 
B,  Fig.  5,  as  follows: 

(i)  Elastic  range: 

KL  «  0.58,     KH  -  0.45,     K^  -  0.78 

Vx  «  0*268  +  0.12  P,     V2  «  0.43H  ♦  0.19  P 

(ii)  fSX as to-plastic  r 

KL  »  0.64,     KM  «  0.50 

\  -  Z   <%s  *  2  *W«     Kep  -  3%fk  CW 

v  «  o.39  h  ♦  0.11  p 

(iii)  uoget 

KL  »  0.50,     KK  •  0»33t     Kj^  *  0.78 

V  -  0.58  a  ♦  0.12  P 

(iv)  Average  value fl  of  the  dynamic  design  factors  in 
elasto-plastic  and  plastic  range  are  given  as  follows: 
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kl  .  Of 6*  +  Q.gQ 
-  0.57 


%  -  ^'/J  »  ^Qg 


-  0.^2  (56) 


L 


irst  xTial  -  j..  1  . /.-Qperties.  Assuming  the  ratio 
of  tensile  reinforcement  to  concrete  areat  p  ■  0.015. 

Assuming     ■  5 
where  *  Jesi,  n  load  ductility  reduction  factor 

•  Ihictility  factor. 
Assumi  ;  C«t  the  ratio  of  maxlmm  resistance  to  the  peak 
load  be  0#7 

jr  •  CL  B 

m    H 

*  .?5  x  60 
■  45  klp» 
The  plastic  moment  resistance  of  beams  with  t     >n  steel 
only  i^  .iven  by  Equation  4.1C,   »f«  5, 

«*  •  Aa  fdy  d  CI  "  fcr?^)  (57) 

,itutinc  the  vr.lues  of  f^f  Agf  and  f'^  in  -ij  na- 
tion 57 i  yields: 

Kp  «  0.688  d2 
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Assuming  Hp  •  ftpB   »  ft^,  Equation  (55)  gives: 
12  MP 


the  values  of  fUt   L  snd  8  gives: 
d  «  9.5  in 

ing  h  *   .5  in  and  d  -  9.5  in 


For  a  value  of  d  *  9*5  in,  Equation  (57)  yields: 
Hp  -  61.92  kip  ft 

Grose  moment  of  inertia  of  the  section  is  given  byi 


4 


*  1160  in 
Bet  moment  of  inertia  is  given  by: 

k3 


net 


.yl' 


+  np 


(X  -  k)2 


v.hurfc 


kd  is  depth  of  neutral  axis. 

Value  of  k  is  s  .1  to  xl  0.42 
therefore, 

Jnet  "  756  *»* 

Average  moment  of  inertia  of  section  is  given 

I.  -  k   (I.  ♦  I—) 


•5  VAg 

908  in 


net" 


"Reinforced  Concrete  Dsslga  Hand  Book',  &m*rie*m  Concrete 
Institute,  second  edition,  p.  54. 
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By  ^uation  (56) 

(16#$)5  (w> 

>0    :ip/f  t 


V; 


■ha  elastic  dc  ion  Y%  1        xven  by: 

«  0,0655 
a  maximum  allowable  displacement,  y     is  n  by: 

17  m  r      e 

«  0.338  ft. 
gfat  of  the  section  is  aqual  to  — *  '7f:>4?oooT'^ 

or 

2.2  kips 

The  mass  of  the  section  is  given  bt,  : 


2 
=*  0.0690  K^f"  to(^ 

(d)  "irst  irial  -  j^uivalent  oystem  •»  ro^er-pies.  The  equiv- 
alent properties  of  the  system  are  as  follows x 

HQ  «  KL  H   -  1.06  kip/sec 

Ke  -  Kn  n       *  0.039  hi   3g2 

1 1 


*a 


nation  (13)   elves: 
(Hj2 


2 

*Te       2 


«  19.30  kip  ft 
The  r  1  period  of  oscillation  is  | :iven  by  the  i^uation 

14(a)    5 


T  -  2 


Jfc 

; 
6 


«  0*056  SQCG. 

(e)    .  i"-.>w      ■>.'•!   1  -    .o.:1.:-    >v>ne  on 

gj ;ity. 

Equation  14  (a)   gives: 

« .  T/_  .  sussr 


*  1.11 


ca 


5i  ■  s§ 


0.75 


♦  7 i       iX  Bt  value  of  m  is  obtained  as: 


T 


-m  ■  0»7*  sec.:. 


3?he  time  Tor  maximum  d  t,   t,  is: 

t^  «  0.71  (T) 
«  0,044  sees. 
e  work  done  ratio,  Cw  is  gives  by  fi  .  6,  ..ndlx  B  as: 

Cw  «  0.32      ' 
Wa$  the  wdwwi  work  done  on  the  equivalent  system  is: 
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.  6.16  ft  kips 

fhe  enerpy  absorbed  by  the  equivalent  system  is  rjiven  by 
ition  6.18f     ieference  5» 

8  "  *me  <Ym  "  °-5  V 
«  7.8  ft  kips 

oince,  £  is  0rsater  than  tff  the  seleoted  proportions  are 
satisfactory  as  a  preliminary  design.*^ 

(f)  xreli^:in>-.:?.y  jjjlfll  £or   i0Rd  Stress.   At  the  fixed-end 
of  the  wall  cover  requirements  results  in  a  smaller  value  of 
d  m   3.75  in  than  at  mid  span  d  *  9*75  in. 
the  fixed  end; 

£hc  estimated  maximum  dynamic  reaction  is  given  by 

V«ax  -  0.5  \ 

*  22,5  kips 

The  value  of  the  allowable  bond  stress ,  ut  is  given  by 
u  -  0.15  **c 

-  450  psi 
fhe  bond  stress %   u  is  given  by 

where 

V  is  the  shear  force  at  the  section  considered,  £  o  is  the 
total  perimeter  of  the  ste^l  bars  and  J  is  the  lever  arm* 
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nation  (58)  yields: 
<     V 

32500 

«  6.6  in, 
AB  •  pbd 

«  0,15  (12)(8.?5) 

-  1.68  in.2 

Trying  iO«  8  round  bars  at  5  in*   spacing  gives 

Ag  »  1.90  in.2  and   £o  «  ?.5  in. 

rherefore,  p  ■  J[ 

=*  0.0182 

At  pinned  end. 

xhe  estimated  maximum  dynamic  reaction,  V^a^*  is  given  by 


1 
■ 
max 


V--  -  V,  H, 


*  13.0  kips 

nation  (58)  yields: 
£  o  *  4.0  in. 
As  «  pbd 

«  0.15(12X9.5) 
-  1.82  in.2 

Trying  No*  8  round  bars  at  5  in*  spacing  gives 
A0  ■  1.90  in.2  and  £  o  *  7.5  in. 
or      p  •  0.0167 
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(    ;    v^ _>■ ■■■   :•   o.'      flttBH  Bid  (:.   -•  UM  l    •>"•-•  ---^c ^_^: 

tion   uy   ..ui^erical   Inte.-.r.-.ti-jn  Method. 

Jubsii   uting  the  values  of  A_f  bt  and  d  in  -<<}uation  (57) 
yields: 

and 

Mpa  -  62.4  kip  ft. 

Ket  moment  of  inertia  is  given  by 


It  .  bd- 


k5^  +  np  (1  -  k)2 


-  12  x  (9.5)5  ^4^~  ♦  0.167  (1  -  0.43)2 
70  in.4 


Therefore 


0.5  (Ig  ♦  \)   •  1015  in. 


nations  (54)  and  (55)  yi«ld  the  dynamic  design  factors 
as  follows i 

iStic  range: 
H, 


' 


8  12s 

E 

30.2  kips 

I? 
870  kip/ft 


i 


0.054?  ft 
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3to-pl  .n tic  ranges 
ep 


- 


5  Ir 

$60  kip/ft 


yep  -  *.  +  Ss-L^ls 

ep 
«  0.084  ft. 


Plastic  range? 


•48.0  kips 


5  *, 


§ 
o 

c 

-p 


•a 


-e  *E  yep 


Deiieeuion,  y(t) 


FIG*  22.  COS  FOR  A  10-3/*"  SLAB. 

Aa  plastic  deformation  ia  permitted 9  the  plastic  Kg  is  de- 
termined by  limiting  the  resistance  to  the  computed  maximum  val- 
ue,  8.  and  equating  the  areas  under  resistance  deflection  curve, 
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•a  under  OABD  -  4  x  3 0.2  x  .0347  ♦  (30.2_>  48) 

2 

x  (0*08%  -  0.0347) 

-  2.435 

Area  OGBD       «  ^  x  48  x  OC  +  (0,084  -  OC)  x  48 
Therefore,   OC  »  I£  •  0*0662  ft 

H  "  *J  -  3^2  •  7a  kip/ft 

3y  Equation  14  (a)  T  is  e^ual  to 


2      I  ...  .   .^;, 


■/ 
■2* 


•  .0^44  sec. 
The  maximum  allowable  deflection,  ▼  ,  is  -iven  byt 

ym  *   ^  P  ye 

«  5  (0.0662) 
»  0.351  ft. 

The  basic  equation  for  numerical  integration  is  niven  by 
Equation  (45) t 

Vi  *yn  (At)2  *  2yn-yn-i 


»]  ere 


in  (At)2  .  c?  -  y(^t)2 

T 
Aasur        ;1m  interval  m  -r— 


0.005  sees. 


yv 


2       (Pn  -  V  25  x  10"6 
Y     (At)     a     * — 6.78  x  0.069 

XI 

«  4.65(lO""4)(Pn  -  R  )   ft.   elastic  range 

2         (Pn  "  V  25  x  10~6 
Yn  (At)  '     " 6.%  COS91 

■  4.65  (10)"^  (P     -  R  )  ft.   for  elasto-plastic 

r  range 

2       (Pn  *  V   (25  x  10     } 
Yn  (  At)t"  "    Kfe  (0.069) — 

■  5*5  10""4  (P     -  S  )  ft.   for  plastic  range 


Table  4. 
reactions. 

Determination 

of  maximum  deflection  and 

dynamic 

t 

sees. 

Pn 
kips 

kips 

P  -R 
n  n 

kips 

V  t)2 

ft 

^n 
ft 

Vl 

kips 

V2 
kips 

0 

60.0 

0 

30.0 

0.01359 

0 

7.20 

11.40 

0.005 

55.0 

11.8 

43.2 

0.02020 

0.01359 

9.66 

15.60 

0.010 

50.4 

34.3 

16.1 

0.00750 

0.04738 

18.90 

18.90 

0.015 

47-5 

48 

-0.5 

-0.00023 

0.08867 

23.70 

23.70 

0.020 

40.6 

48 

-7.4 

-0.00406 

0.12973 

23.00 

23.00 

0.025 

35.9 

48 

-ll.l 

-0.00610 

0.16673 

22.40 

22.40 

0.030 

31.0 

48 

-17.0 

-0.00935 

0*19763 

21.82 

21.82 

0.035 

27.2 

48 

-20.8 

-0.01142 

0.21918 

21.36 

21.36 

0.040 

21.4 

48 

-26.6 

-0.01460 

0.22931 

20.67 

20.67 

0.045 

16.6 

48 

-31.4 

-0.01722 

0.21475 

20.0 

20.0 

0.050 

11.7 

48 

-36.3 

-0.02010 

0.18297 

19.50 

19.50 

<Vma*  "  °'22931  -  0.23  ft. 
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;  >-  c  loulatieaa  for  I  ili  %i 

for  t  •  0.025  and  n  «  p,  Ir  is  ^iven  b,  r 
T5  .  T4  (At)2  +  2  Y4  -  X? 

-  -(0.00406)  ♦  2  (0.12973)  -  0.08867 
«  0.16673  ft. 

Hince  Yc  «  0.16673  >  0.0662,  R^  «  48 

The  dynri&ic  reaction  is  given  by  Equation  (55)  and  is  e  iual  to 
V  •  0.38  H  ♦  0.12  P 

-  0.38  (48)  +  0.12  (35.9)  -  22.40  kip 
By  Table  4f 

Maximum  deflection  »  0*2293  ***  which  is  less  than  allowable 

(h)   w  r  and  Bond  Strength*  By  Table  4f  Vaax   *  23.70  kip 
(i)  Tor  no  shear  reinforcement,  allowable  shear  strength  is 
given  by 

? l       0.04  f »c  +  5000  p 

•  0.04  (3000)  ♦  5000  (0.0167) 
-  120  +  83.5  *  ^03.5  P«i 


Shear  stress  v  « 


m  •  %&l?.hj  -  258*° p8i 


As  the  shear  stress  is  sore  than  allowable,  shear  rein- 
forcement is  required  for  (258  -  203«5)  •  54.5  psl. 

Contribution  of  shear  reinforcement  to  allowable  shear 
stress  is  civori  by: 
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f 


2U5. 


r  «  -?r-  *  ,^\<,   =  0*00156 
I 


drying  one  SF&«  3  round  bar,  Afl  ■  J.  11  in. 

s  *  8.10  in.,  using  /3  bars  at  8  in.  spaclr- 
At  the  top  of  the  wall: 

V  -  203*50  psi 

X      25700 

«  243  pel 

Using  i o.  3  round  bar  at  S  in.  c/c 

Summary!  Slafc  10.75  in.  thick,  p  -  0.0163  she  tnt 

for  bottom  5  ft.  and  top  5  ft.  1  So«  3  *•  3  in.  spacing. 

Design  of  Soef  Slab 

The  roof  slab  Is  designed  a*  a  one-way  reinforced  concrete 
■lab«  The  slab  is  permittf     deform  into  plastic  range  devel- 
oping plastic  hinges  at  the  support.,,   £na  critical  roo  "-.  lab 
loading  is  the  incident  overpressure  versus  tine  curve  shown  in 
the  Fie".  lfc«  This  loading  results  fro®  the  1     rfave  moving 
parallel  to  the  long  axes  of  the  building. 

The  blast  loading  curve  of  the  Fig.  16  is  i    .zed   to  a 
triangular  loud  of  10  ..  ;e  and  the  duration  of  0*38 

sacs,   following  the  dynamic  design  procedure,  the  thickness  of 
slab  was  determined  to  be  7.25  in.  with  p  m   0*016*   lauaerical 
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integration  met'-  od  was  used  to  check  the  design  and  to  find  the 
dynain  ic  re  *  Q  t  i  0  D  s  • 

•eliminary  Column  Oesi($n 

In  deten        •  spring  constant,  the  column  height  is 

taken  e^ual  t-   .75  ft.  from  center  line  of  girder  tc  the  top 

oi  the  xootii  .     \   resistance  computations  are  based  on  the 
cle. r  height  of  1$   feet. 

(a)  Design  loading*   In  computin ;  the  total  eenoentx 
load  on  the  frame f  it  is  assumed  that  the  wall  slab  transmit 
the  blast  e  ,uaily  to  the  roof  slab  and  the  foundations. 

3?he  design  load  as  idealized  froza  the  computed  loading  as 
shown  in         g  1?  is  defined  by: 


/-~\ 


x1   =     .-  '- 


v  •  (sees.) 


H 


FIG,  2$.      XBEALIZS9  LOAD  3tt 

(25.5:  (1^(16)  (16/2) 
1C 

466  kips 

0.062  »«0  . 

m 


14.6  kip  i;ec 
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(b)  Mass  OoiSi)Utation» 

Weight  of   the  roof  slab  ■  (2a^Sj|JiQ.)  (^6  ffljjjl*'*) 

«  43.8  ki  t? 

Weight  of  the  girder  (assumed)   «  ^jfff  j&   .^^2) 

»  19 • 2  kips 
Weight  of  the  three  columns  (assumed)   »  i;  •  ■     .-i1^q^^  ^ 

■11*7  kip© 

12   UOOOX 

-  57.0  kip 

10 
The  mass,  m,  of  the  sin  le  degree  freedom  system  is  riven  by 

m,   Is  e  gusJL  to  th-  j        ss  of  (roof  ♦  fch<  tier)   ♦  1/5  mass 

of  the  (columns  plus  wall  slabs) « 

there.-*  W 

*"»  « ~»  fci.P.  s  e  c 

m   «    2 •  S3   "f* Mir"' TV-  *"  — 

(c)  First  vri&l  ctual  Properties* 

-naming  p  •  0.015, 

•  6, 
and  GR  «  0.50 

\  m   °R  B 

*  233  kips 

The  maximum  design  moment,  fU   ,  under  axil  load  is  given 
by  equation  (7*15)  Reference  3 
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«D  "  -S-H  (59) 

.  9 

where 

n  •  number  of  columns 

nation  (59)  yields: 
Up  «  5*30  kip  ft. 

Average  roof  overpressure,  frc  <  i  .  20  *  0.5  x^si. 

je  blast  load  per  column,  due  to  roof  overpr     e  is 
equal  to 

3  C1C0Q) 
-  16?  kips. 

.  1  per  column  is  given  by 
1/3  (43«8  ♦  19.2) 
«  23  :  i 

Column  design  load  P^  «  167  +  23 

*  190  kip* 
2he  moment  for  failure  in  a  symmetrically  reinforced  column 
is  given  by  the  Equation  (4*32)  reference  9* 


MD  "  As  %  *'  +  *B  <f  "  i.?  b^'dc  <«» 


id  c;e 


Pn  •  load  on  the  column 

t  «  thickness  of  the  column 
d*  •  distance  between  centrolds  of  the  compression  and 
the  ten  .  ile  steel 
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Assuming  b  ■  12  in. ,  p  •  Q.G1  and  substituting  the  valuea 
of  iDf  KD,  fdc$  and  f^  in  the  Equation  (60)  t  and  solving  for 
t  yields: 

t  «  25.5  in* 
Trying  t  •  24  in. 
Equation  (60),  for  t  *  24  in.  yields: 

ftp  »  516  kip  ft. 


~: 


i 


T 


-j 


i_   ± 


J 


«P 


FIG,  24.  SEGVXQM   OP  A  GOLUHH 


For  |"«  §J^  -  0.105, 

and  p  «  p»  »  0.015 

k  is  given  ■  e  ual  to    0.30 


by: 


#et  moment  of  inertia  of  a  doubly  reinforced  beam  is 


2 

♦  np  (1  -  k)~|   (61) 


I.    •  fed-7 

r,      dfT 

Equation  (61)  yields: 

It  m  1090 

0  in,4 

Gross  moment  of  inertia  is  given  by 
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Xg  "if" 

-  13824  in.4 

Average  aoaent  of  inertia  Ifl  obtained  as 

Ia  -  3     in.4 

£h*  spring  constant  of  the  colon  la  C?»3; 

iiefsresce  1 

aation  (62)  yields: 
k  -  2880  kip/ft. 
For  H3  •  i>16  ft      L.  ,      [OfttlOB  (59)  yields: 
^  •  2$8  kip/ft, 

the  *;ic  def  onf   X     1  vein  by, 

*e  -r 

18  •  ■...■•   rti 
and  the  aaxiauc  allowable  deflection  is  given  bg  I 

■  .498  ft* 
natural  frequency  is  given  byi 

t  •  2 

m   j.197  sees. 


(d)  Jfirst   ^rial  -    ;o?k    >one   Versus  J^ner?ar  absorption 

i*j cit^ . 

i'he  ratios  ^-  and  Ca  are  given  by 


n 


C   o7l5? 

n 

and  Cjj  * 

-  0.51 

% 

v   this  ratio,  value  of  ^p  is  obtained  frm   I  ..  7f  Appen- 
dix B|  and  is  equal  to 

ig*  •  1.04 

or 

t  »  ...0645  sec. 

The  original  assumed  load-tiiue  curve  shown  in  Fig.   23  is 
revised  to  obtain  a  closer  appro:         Aon  up  to  the  tii^e  t   . 
The  impulse  up  to  t  -  0.10  is  H  •  1.167   (from  Pig.   1?.) 
Therefore 

T  -  2| 

»  0.0923  sec. 


T     . 


u.47 
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For  S-  »  0.47  and  CR  -  0.51» 

gP-  1.023  it*om     &i  .  ?»  ppendix  B) 

or 

tm  «  0.095     . *  wMch  is  O.K. 

The  work  done  ratio a  Cyf  is  obtained  from  the  ftg«  6, 
Appendix  B  and  is  equal  to 

V  is  given  by 

VI  .£ 

S  ■  aS 

*  86 »0  ft  kipe 

9  energy  absorbed  by  tl.         ;nt  system  is  £.;iven  by 
the  Equation  6.18»  Reference  5  as  follows s 

S  -  Hm  (xm  -  °*5  xe) 

*  107  ft  :-i.  3 

£   W%  hence     n  is  satisfactory  as  a  preliminary  design. 
2ho   design  was  also  found  to  be  adequate  after  checking  by 
the  numerical  integration  nethod. 

Design  of  itoof  Girder 

The  fre&e  of  this  building  consists  of  three  rectangular 
columns  supporting  a  rectangular  girder  which  forms  a  tee  bean 
;h  thfl        lb«      ^oof  girder  i       nod  to  resist  the 
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FIG.  25.   SECTIONS  OF  THE  GIRDER  AT  THE  MID-SPAN  AND 
AT  THE  SUPPORT. 
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combined  vertical  loads  on  the  roof  and  lateral  loads  on  the 
frame.   Ithough  the  other  structural  elements  of  this  building 
are  permitted  to  deflect  plastically,  the  girder  is  deuianed  to 
behave  elastically  so  that  pro;>er  restraint  is   maintained  for 

the  column  throughout  the  deflection  of  the  fr  e. 

11 
1'he  design  bending  moment  of  the  girder  is  the  sum   of 

moments  due  -co:  (a)  the  roof  slab  dynamic  reactions,  (b)  the 

static  loads,  and  (c)  the  frame  action.   The  moment  due  to 

frame  action  is  e ^ual  to  one-half  the  column  plastic  moment. 

The  girder  is  designed  as  a  tee-beam  in  the  ic  ion  of  ps- 
Itive  3  0iic,  .  in  the  region  of  negative  moment  near  the  inte- 
rior support,  the  girder  is  designed  as  reet  ion. 

The  sections  (  I  •  25)  are  found  to  be  adequate  for  this 
buildin ;  after  checking  them  by  the  numerical  integration  meth- 
od. 

C0KCLU3X0BS 

Pictures  designed  to  resist  the  blast  loads  are  subjected 
to  completely   LfJ    b  types  oi  loads  than  those  considered  in 
conventional  design!  .    ie  to  the  large  magnitude  ar    MBit 

racter  of  loading,  the  di  =j  are  based  on  dynamic  analysis. 
In  order  to  simplify  the  analysis,  a  given  structure  is  replaced 
by  a  dynamically  equivalent  system. 

e  design  example  presented  in  this  study  is  analyzed  by 
the  I       1  integration  method.    j  to   economical  reasons, 
the  aembere  c        ticture  are  allowed  to  deflect  plastically 


b€ 


and  the  dynamic  yield  stresses  are  used  for  design.   It  is  felt 
that  concrete  members  of  greater  mass  are  more  suitable  for 
blast  resistant  construction  than  steel  members  of  smaller  mass* 
es  due  to  the  inertia  of  members • 
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FIG.  6.   WORK-DONE  RATIOS  FOR  TRIANGULAR  LOADS  AND  ONE- 
DEGREE  SYSTEM.   (P.  53,  Ref.  8.) 
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APPENDIX  C  -  SYMBOLS 


A    Area  of  tension  steel  in  reinforced  concrete  menber. 
s 

B   Peak  value  of  externally  applied  load. 

R_ 
Cg   ^atio  of  maximum  resistance  to  peak  load,  Gfi  ■  -§c 

Cm   Ratio  of  load  duration  to  natural  period  of  oscilla- 

T 
tion,  Cm  -  -m- 

n 

C,,   Ratio  of  maximum  work  done  to  absolute  maximum  work 

W 
done,  Cw  «  7T- 

P 

X 

DLP   Dynamic  load  factor  *  y" 

a 

E   Energy  absorbed  by  the  equivalent  system  f  modulus  of 

elasticity. 

f.    Dynamic  yield  strength  of  steel  (psi). 

f • d    Dynamic  ultimate  compressive  strength  of  concrete. 

H   Impulse  per  unit  area,  impulse  per  foot  of  length. 

H    Equivalent  impulse  acting  on  equivalent  system. 
e 

I    Average  of  the  gross  and  transformed  moment  of  inertia. 

I    Moment  of  inertia  of  gross  section. 
S 


Moment  of  inertia  of  transformed  section. 


K~   Load  factor 


KrM       Load  mass  factor. 


KM 


Mass  factor. 


77 


K«  assistance  factor. 

k  Spring  constant* 

k«  ISffectivs  spring  constant* 

ke  Equivalent  spring  const  rant. 

k*»~>  Spring  constant  in  elaato  plastic  rang** 

KD  Haxiaum  design  moaent  in  a  member  under  axil  load  »,. 

Mp  Plastic  bending  moment  under  bending. 

M~  I  lastic  bending  moment  at  center  line  of  beam  or  slab. 

Ila-  Plastic  bending  moment  at  support. 

p   P 

m  Haas  per  unit  length  (kip  -  see  /ft  )• 

..  Has.  «   e,uivalent  8y£te»  CM.  -  -A**>. 

1  *  Haximum  axil  load  on  column* 

Pa  Overpressure  existing  in  incident  shock  wave  for  any 
value  of  t  -  t  .,• 

P  Initial  peak  incident  overpressure* 

p(t)  Actual  load  on  a  structural  element  as  a  function  of 
time* 

p  Katio  of  tensile  reinforcement  to  concrete  area  » 

3  Total  resistance  of  structural  element* 

&m  Maximum  resistance  developed  by  a  structural  system* 

B^  Maximum  resist  nee  in  equivalent  system* 

T  Duration  of  external  loads* 

Tn  Katural  period  of  oscillation. 


?8 


t         Timo  required  for  maximum  displace/. ent  of  element. 

t    Duration  of  the  positive  phase  of  incident  shock  wave. 
o 

t     ,;e  interval  used  in  numerical  analysis. 

7   Dynamic  reaction. 

v    near  stress. 

Maximum  work  done  on  equivalent  system  by  e^uivilent 

load. 

Fictitious  maximum  work  done  on  equivalent  system. 


x    Maximum  displacement. 

yft   Reflection  of  equivalent  system,  limiting  elastic 

deflection. 
y     Limiting  deflection  in  elasto  plastic  range. 

Maximum  displacement. 

yn   Displacement  at  time,  tn. 

<    Design  load  ductility  reduction  factor. 
p       Ductility  ratio. 
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